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result that the essential dynamic part of the Hamiltonian for the system of identical particles is 
the Reduced Hamiltonian operator describing relative movement of two particles inside the sys- 
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I. INTRODUCTION 



The wave function of the self-bound system in the absence of external fields must be in- 
variant in respect of translations in space. This invariance of wave function means that it is 
dependent only on intrinsic degrees of freedom of the system. However, the traditional meth- 
ods for quantum systems description, such as Shell Model or Hartree-Fock Self-Consistent 
Field (SCF) method, produce wave functions dependent on a set of one-particle variables, 
hence also on the system center of mass radius- vector. This shortage of mentioned methods 
is well-known, however, the wave functions dependent on one-particle variables are very 
attractive because they allow simple procedure of antisymmetrization. In some cases this 
approximation with the Hamiltonian and the wave functions, both dependent on redundant 
variables, cannot produce any serious problems. The moving and exciting in uncontrollable 
way center of mass of the system, such as an atom, molecule, or electron gas in solid state, 
is not the problem for the system under consideration (i.e., the system of electrons) but the 
problem of external fields, keeping these electrons together to form such a system. However, 
for nuclei and hadrons the translational invariance appears to be a real problem. These 
systems are essentially self-bound. The external fields, able to dictate behavior of such a 
system in normal conditions, arc very weak in comparison with interactions of nucleons or 
quarks, i.e., the structural elements of nuclei or hadrons. This invariance of corresponding 
wave functions creates less investigated additional correlations between particles. 

A quantum system without correlations is a system, whose wave function can be present 
as a product of wave functions of subsystems, composing the entire system. Such factor- 
ization of wave function is possible when interactions among these subsystems are absent 
or can be described by some central field forcing these subsystems to compose the system. 
Such systems do not exist, this is only a more or less acceptable model simplifying the 
description. If particles compose a system, they always are correlated, and first of all due 
to interactions among them. However, besides dynamic correlations, there exist additional 
correlations, characteristic namely of quantum systems. These are correlations created by 
different symmetries of the system. Any symmetry has a corresponding operator commuting 
with Hamiltonian. Wave function of the quantum system must be an eigenfunction of these 
operators. This introduces the new quantum numbers characterizing the system, but at the 
same time requires linear combinations of mentioned products of subsystem functions, hence 
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additional correlations. 

The first among them is invariance of the system in respect of rotations and inversion. 
The corresponding quantum numbers are angular momentum and isospin, projections of 
these quantum numbers, and parity. The second invariance is invariance of the system in 
respect of identical particle permutations. For the system composed of identical fermions, 
this symmetry is known as the Pauli principle. The latter kind of correlations is created by 
the mentioned above translational invariance. At any kind of correlations taken into account, 
one gains some quantum numbers for the entire system, but at the same time particles lose 
individuality and no one of them can be described by a complete set of necessary quantum 
numbers. For example, binding the momenta one loses momentum projections for individual 
particles. Antisymmetrizing the function, the determinants appear, states for individual 
particles are not defined at all. Only configuration, i.e., a set of restricted quantum numbers 
characterizing the entire system, is defined. For translationally invariant function one cannot 
define even the configuration. 

Some recent methods, developed for few-body systems with strong interaction, are de- 
signed to produce translationally invariant wave functions. These are methods based on the 
Faddeev or Faddeev-Jakubovsky equations on the Green Function Monte Carlo (GFMC) 
method or on the expansions in a large basis of harmonic oscillator functions [3i] . In the 
first case, the translational invariance of wave function arises from equations written in 
intrinsic variables. For GFMC the translational invariance is given by infrastructure of vari- 
ational wave function, spatial part of which is composed of Jastrow multipliers dependent 
on translationally invariant differences of one-particle coordinates. The mentioned latter 
calculations operate with the harmonic oscillator wave functions as basis for realistic cal- 
culations. These functions, if complete basis of ones corresponding to a given number of 
oscillator quanta is taken into account, can be projected to superpositions with unexcited 
center of mass of the system. 

The goal of the current study is new efficient method of the construction of the trans- 
lationally invariant wave functions for the system of identical fermions with arbitrary po- 
tential of pair interaction. These functions are eigenfunctions of intrinsic Hamiltonian, i.e., 
Hamiltonian dependent on intrinsic variables. One ensures translational invariance for wave 
functions taking as arguments the set of translationally invariant spatial variables, so called 
Jacobian variables, defined as in jj] . The formalism is free of any phenomenological param- 
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eters (as variational parameters in GFMC) or fields (as in any model exploiting the shell 
model idea) as well as of any "effective interactions". 

The basic point of the method is the well-known result that the essential dynamic part 
of Hamiltonian for the system of identical particles is the Reduced Hamiltonian (RH) oper- 
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6|. The system of eigenfunctions of this operator takes into account all correlations 
caused by interaction. The Schrodinger equation for RH operator is very simple, not more 
complex than the Schrodinger equation for relative movement of two particles interacting 
by given potential. This really looks like a one-particle problem, like in the Shell Model. 
These eigenfunctions are basic in our formalism to the wave function construction. The 
main problem in this approach appears to be antisymmetrization of a component (i.e., the 
construction with lower degree of antisymmetry) of wave function. The developed univer- 
sal algorithm for antisymmetrization gives us a possibility to perform this operation in a 
simple way and keep numerical approximations under control. In the mentioned basis with 
correlations taken into account in advance, convergence is fast. 

In Sec. 2 we present the basic definitions and introduce the problem under consideration. 
Sec. 3 is devoted to description of the new original procedure for antisymmetrization of 
wave functions for translationally invariant systems composed of identical fermions. In Sec. 
4 one considers the entire problem based on factorization of the Hamiltonian matrix. The 
summary of the developed method, numerical test and conclusions are given in Sec. 5. 

II. CORRELATED COMPONENTS OF WAVE FUNCTION 

The suggested consideration is based on the presentation of an intrinsic Hamiltonian for 
the system of N identical particles as a sum of A^(A^ — l)/2 two-particle operators: 



^2 ^ 



1=1 ij = l{i<j) 



N 



i,j=l(i<j) 



N 



E -^(V.-V,)' + y(r,-r„a,a,) ^ J] (1) 



«j=i{«<i) 



where V {ri — rj,aiaj) is the potential of interaction, dependent, as usual, on difference of 
radius- vectors of interacting particles i and j and on discrete degrees of freedom cTj, aj of 
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these particles, such as spins, isospins, etc. The term fv^Vl /2mN is the kinetic energy 
operator of the center of mass of the system. 

The arbitrary eigenvalue of this Hamiltonian, due to antisymmetry of corresponding 
eigenfunction, equals the number of terms present in the sum, multiplied by the expectation 
value of one of these operators, i.e.. 




N 



(h^j) = {H,j) . (2) 



, N . ^ 

Here J j — N {N — 1) /2 is a binomial coefficient and ^ij — — rj) is one of the 

Jacobian variables. To simplify the expressions, let us take defined values for indices, namely 
i — N — 1 and j — N and mark the corresponding Jacobian variable as ^ = ^jv-i n- 

The 

operator Hn-i,n = H is called the Reduced Hamiltonian operator and is the main in 
further consideration. The best way to eliminate dimensions of physical variables in above 
expressions is to present the multiplier h^/m in the form 

— ^nujxh\ (3) 
m 

where hiij can be considered as a parameter of energy and 6 as a length parameter. Really, 
one free parameter appears here, uj or h. Let us consider below the dimension-free values for 
energy £ = E/fiu and variable ^ = ^/h. The value for parameter h without any problems 
can be taken equal to one. Finally, the dimension-free Reduced Hamiltonian equals 




(4) 



Thus, the expression for the eigenvalue of Hamiltonian in terms of hjj is £ — (H) . 

This is a very interesting result stating that calculation of eigenvalue does not require the 
total Hamilton operator. Only the significant part of the many-particle Hamiltonian with 
pair interaction - the two-particle RH operator - is necessary. However, practical application 
of simplifications caused by this observation is not easy. Due to noncommutation of H and 
7i operators and antisymmetry of wave function one cannot prescribe some set of quantum 
numbers for the state of the last two particles. The way of solving this problem was found 
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years asro, when coefficients of fractional parentage (CFP) for the atomic shell model were 

un 

introduced [3, |8|. These are the coefficients of the antisymmetric wave function expansion in 
terms of linear combinations of some easier structures with a lower degree of antisymmetry, 
later on called components of wave function. For the problem under consideration the 
components are products of antisymmetric function for the last two particles dependent 
on variables present in the given Reduced Hamiltonian operator, and the antisymmetric 
function dependent on a set of all residual variables ^ . In other words, the idea of coefficients 
of fractional parentage realizes the antisymmetrization in space of quantum numbers of the 
complete system of components. Let us mark the antisymmetric wave function for the N 
particle system as 

^^AA/(l,2,...,iV-l,iV), (5) 

where S is the corresponding eigenvalue of Hamiltonian, A is a set of "good" quantum num- 
bers, such as momentum, parity, isospin, etc., and M is a set of corresponding projections 
of momenta quantum numbers. 

The mentioned construction with a lower degree of antisymmetry, i.e., the component of 
wave function, is: 



*(rA,.p)AM (1. 2, iV - 1, iV) = ^ v&f am (1, 2, iV - 2) {N - 1, N) 

(6) 

It is a bound momenta function with exact quantum numbers AM coinciding with ones for 
complete wave function. The function \E' present on the right-hand side of equation is a 
"spectator" function dependent on variables of ffist {N — 2) particles. AM is a set of corre- 
sponding quantum numbers and its projections. F marks a set of all the remaining quantum 
numbers necessary to ensure completeness and orthonormalization of basis of "spectator" 
functions. The second function is eigenfunction of the RH operator, hence the best way 
is to choose the quantum numbers p and /i as a complete set of eigenvalues of operators 
commuting with the RH operator, e marks an eigenvalue of the Schrodinger equation for 
RH. The last factor is the product of Clebsh-Gordan coefficients for momenta and Kronecker 
deltas — for parities, isospin projections, etc. 



A p A 

M n M 
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The most interesting for us are coefficients of the wave function expression in terms of a 
complete set of components: 

^eAMil,2,...,N-l,N) = Y,{£A\\TA,ep)^^^j^^^^^^,,il,2,..-N-l,N). (7) 

TA,ep 

The system of equations for these coefficients consists of the Schrodinger equation and 
condition that wave function of the system, composed of identical fermions, must be an 
eigenfunction of the antisymmetrization operator corresponding to the eigenvalue equal to 
one: 

{n - 8) ^^AAf (1, 2, iV - 1, iV) = 0, (8) 

<^sKM (1, 2, iV - 1, N) = 0. (9) 
Obviously, the antisymmetrizer, present in Eq. Q, must be properly normalized, i.e.: 

^=jr^ E (10) 

where P are permutation operators of symmetric group Sn and 5p is the parity of permu- 
tation P. Namely due to this normalization the antisymmetrizer is a projection operator, 
i.e., satisfies the condition AA=A. The same condition for matrix of this operator simplifies 
further consideration. 

However, simultaneous solution of both equations is very problematic. As usual in quan- 
tum mechanics, ffist of all one needs to construct a complete set of solutions for a simpler 
equation and later on, applying these functions, to construct solutions of a more complex 
equation. The simpler equation in our case is the second equation, so ffist of all it is necessary 
to antisymmetrize the components. The antisymmetrized basic functions, eigenfunctions of 
antisymmetrizer, are: 

vl/rAM (1, 2, ...,N-1,N)=J2 (rA||f A, ep) $(rA,.p)AM (1' 2, iV - 1, N) . (11) 

rA,ep 

Namely coefficients, present in this expansion, are the coefficients of fractional parentage; 
r is an additional quantum number necessary to mark the antisymmetric basic functions. 
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The last operation is diagonalization of the Hamiltonian operator while applying this basis 
of functions. For this one needs to define the expansion of wave function 

*£AM (1, 2,...,N-1,N)^J2{S\T) *rAM (1, 2, ...,N- 1, N) (12) 

r 

and to solve an algebraic eigenvalue equation for these coefficients: 

J2 [(r \n\ V) -s]{s\ V) = o. (13) 

r' 

The final expression for wave function in terms of components with a lower degree of 
antisymmetry looks as: 

(1, 2, TV - 1, TV) = ^ I r) J] (rA||f A, ep) $(rA,ep)AM (1> 2, N-1,N). 

r rA,ep 

(14) 

III. ANTISYMMETRIZATION IN CORRELATED BASIS 

Let us first of all formulate the well-known results, concerning antisymmetrization of dif- 
ferent wave functions in a way giving possibility to generalize it for translationally invariant 
functions. 

The best known procedure is antisymmetrization of many-particle function dependent on 
one-particle spatial variables and intrinsic degrees of freedom, when particles do not interact, 
correlations are absent, and wave functions split into product of one-particle functions. In 
this case, antisymmetrization can be performed in given configuration, i.e., for each set of 
one-particle states provided that coinciding one-particle states are absent. This construction 
equals the determinant composed of mentioned one-particle functions with permuted one- 
particle variables. 

The simplest fractional parentage expansion for this wave function is expansion for de- 
terminant in minors, i.e., in constructions antisymmetric only in respect of permutations 
of the first {N — 1) particles. The last particle is not antisymmetrized with particles com- 
posing this "spectator" system. However, a state for this particle in an antisymmetrized 
function runs over all states present in configuration under investigation. The coefficients of 
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expansion for antisymmetric function, equal to the determinant of the A^-th order in terms 
of minors of (A^ — 1) order, are the one-particle CFP. They are well-known from linear alge- 
bra. Due to normalization of wave function and all minors they equal ±l/v^. In this case, 
definition for the two-particle CFP is completely analogous. They are coefficients present in 
the same determinant expansion in terms of minors of the (A^ — 2)-th order multiplied by 
antisymmetric wave function for remaining two particles. 

The antisymmetrization and fractional parentage expansion for more complex function 
with bound momentum-like, i.e., orbital, spin, isospin, etc., quantum numbers, can be for- 
mulated in analogous way. In this case, there appear shells, i.e., groups of particles having 
the same one-particle radial wave function but different projections of orbital and intrinsic 
momenta of particles. However, the antisymmetrization for such a function as in previous 
case is impossible. The operation of momenta binding results in sums over projections of 
quantum numbers, hence a complete set of one-particle variables for every particle stays 
undefined and the old scheme of antisymmetrization cannot work. Namely the antisym- 
metrization in the defined shell creates traditional coefficients of fractional parentage. As 
in the previous case, they are independent of one-particle radial wave functions, hence are 
universal. For more complex systems, containing more than one shell, antisymmetrization 
is not so complicated as believed because radial functions for different shells are orthogonal, 
and this significantly simplifies consideration. The CFP, defined for such a complex config- 
uration, are expressible in terms of traditional CFP and transformation matrices. Again in 
this case antisymmetrization can be performed for defined configuration, it is in finite space 
of components basis. 

Antisymmetrization for translationally invariant wave function is a very complex problem 
in comparison with considered ones. In the form given below, it is introduced in and 
still applied only to three-dimensional harmonic oscillator functions This is the simplest 
application of the proposed scheme because only for harmonic oscillator basis this operation 
can be performed in finite space of basic component functions. The reason for this is a 
very simple structure of the many-particle harmonic oscillator Hamiltonian, invariant in 
respect of any permutations of Jacobian variables. This symmetry creates a conserving 
quantum number - a total number of oscillator quanta proportional to the eigenvalue of 
such a Hamiltonian, and the possibility to perform antisymmetrization separately for each 
value of this integer number. However, this simplification is not so useful as it looks at first 
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glance because in many cases long expansions in this basis are necessary and convergence of 
expansions is very slow. 

The procedure one applies is defined in a universal way and is based on calculation of 
the antisymmetrizer matrix in components basis and spectral decomposition of this matrix 
in case when the eigenfunctions of the RH operator are applied as "actor" functions. The 
formalism is devoted to bound states of the quantum system description, so in expansion 
of component one can apply square- integrable functions for "spectator", as well as square- 
integrable functions for "actor". Obviously, in most interesting applications RH has both 
discrete and continuous eigenvalues. The necessary basis for the corresponding Sturm- 
Liouville problem can be constructed taking proper boundary conditions for eigenfunctions, 
for example, taking all eigenfunctions equal to zero at some value of argument larger in 
comparison with a characteristic measure for the system dimensions. 

The next problem is a complete system of functions for " spectator" . For these functions 
one cannot formulate any equations, hence one is free in choice. The simplest translation- 
ally invariant basis, as mentioned above, is basis of three-dimensional harmonic oscillator 
functions, enriched by intrinsic degrees of freedom of particles. 

Once the construction of basis for components is completed, one can start calculating 
the matrix of an antisymmetrization operator. If, as defined, one has antisymmetrized 
"spectator" as well as "actor" functions and is going to apply this operator to the matrix 
calculation, the best factorization for the antisymmetrization operator is as follows j^: 

Ai^2,...,N = ^Ar-l,Af^l,2,...,Af-2^1,2,...,Ar-2;Af-l,Ar^l,2,...,Ar-2^Ar-l,Ar, (15) 

where 



yi,2,-,N-2;N-l.N = \ [1 — (A^ — 2) Pn-2,N + Pn -"i.N -iPn ~2,n\ (16) 

2 / \ 2 

and Pjj are transposition operators. 

The normalization of antisymmetrizer given in Eq.f llOl) makes corresponding matrix idem- 
potent or projection matrix: 



AA = A. (17) 
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This matrix is real and symmetric (A+ = A). Obviously, its eigenvalues equal ones and 
zeroes. Hence, the sum of eigenvalues of this matrix equals the number of eigenvalues 
equal one. By definition, it coincides with trace of matrix A and defines the rank r of this 
matrix. Eigenvectors, corresponding to zero eigenvalues, do not have a physical meaning. 
Antisymmetric states can be described only by eigenvectors corresponding to eigenvalue 
equal one. This matrix has an attractive feature - each entry Oj^ of this matrix equals the 
product of the i-th and j-th rows (columns) of this matrix. Hence, each diagonal entry 
of matrix equals the sum of squares of entries of this row or column. This is a useful 
measure of projectiveness of any calculated part of matrix. Moreover, if some diagonal 
entry of this matrix equals zero, all entries of corresponding row and column are identical 
zeroes. This means that such a basic component must be removed from basis because it has 
some additional symmetries and cannot be antisymmetrized. The second conclusion from 
this feature of the projection matrix is that if some diagonal entry of this matrix equals 
one, again all remaining entries of corresponding row and column equal zero. This means 
that the corresponding basic component is already antisymmetric and needs no additional 
antisymmetrization. 

The spectral decomposition for this matrix is 

A = FF+, (18) 

where F is the matrix composed of eigenvectors, corresponding to unit eigenvalues, as 
columns. Every column of this matrix defines antisymmetric wave function for parti- 
cles. Matrix F has n rows {n equals order of matrix A) and r columns. Every column is 
normalized and they are orthogonal, i.e.: 

F+F = 1. (19) 

The entries of matrix F are namely the coefficients of fractional parentage for any of 
mentioned components bases. For this matrix construction, observation that every column 
of matrix A is eigenvector of this matrix, corresponding to the eigenvalue equal one, is very 
useful. However, the number of these eigenvectors exceeds the rank, hence some of them 
are lincarly-dcpcndcnt. Moreover, they are not orthonormal vectors, but direct calculation 
of matrix F is not the problem due to the well-known result that each symmetric matrix of 



11 



an order, equal n, and a rank, equal r, can be written in a form 

A= I ^ I , (20) 

' Q QZ-^Q^ ' 

where Z is nondegenerate submatrix of the order r. The matrix of spectral decomposition 
for matrix A is: 



(21) 




The normalization condition for F looks as 

F+F = Z+Z-1/2Q+Q2-1/2 _ i_ (22) 

It can be written in the form 

Q+Q = Z(1-Z). (23) 

IV. EIGENVALUES OF HAMILTONIAN MATRIX 

Once the matrix F is constructed, one is ready to obtain a matrix of the total Hamiltonian. 
As it was mentioned, first of all one must take into account that this matrix due to identity 
of particles equals the matrix of the RH operator. Having in mind that the coefficients of 
the antisymmetric function expression in terms of components are CFP, i.e., entries of F 
matrix, one can present the matrix of the Hamiltonian in a form 

F+HF. (24) 

Here the matrix of the RH operator is given in components basis. This is our goal because in 
this basis matrix H is diagonal. The best choice is to rearrange the entries of this matrix in 
the nonvanishing order. By the way, this ordering of basic functions looks very natural and 
means that one takes the function with " actor" being in the lowest possible state as the first 
basic function, then in the first excited state - as the second basic function, etc. Obviously, 
for every "actor" function one has to take a necessary set of "spectator" functions. 
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The order of the Hamiltonian matrix equals rank r of matrix of an antisymmetrizer, 
but the order of the diagonal matrix of the RH operator is larger, equal to the order of 
the antisymmetrizer matrix n or, in other words, the number of basic components taken 
into account. This infrastructure of the Hamiltonian matrix is very useful in cases when 
dimensions of matrices present in expression are very large. Taking into account this in- 
frastructure, one can construct different approximations for eigenvalues and eigenf unctions. 
The idea is as follows. The antisymmetric states, involved in the expansion, are given by 
columns of matrix F. Let the rank of this matrix be equal one, this means it contains only 
one column. In this case the eigenvalue of the Hamiltonian equals: 

£ = Eifl^ + £2/21 + ••• + ^nfll = eiUJl + e2UJ2 + ••• + £nUJn- (25) 

Here are the eigenvalues of the RH operator arranged in the nonvanishing order. Due 
to reality of matrix F the squares of entries are the probabilities (due to normalization 
this sum of squares equals one, i.e., Y17=i fii ~ ^) '^^ different states of the RH operator 
in the expansion for the energy eigenvalue of the total Hamiltonian. The expansion of the 
same kind can be obtained for each eigenvalue after matrix diagonalization because the 
corresponding eigenvector equals a linear combination of columns of matrix F: 

n n 

£ = '^eiUJi, ^u;i = l, uji>0. (26) 

i=l 1=1 

As it follows from definitions, all Ui are nonnegative. At the same time some first eigen- 
values of RH are negative and remaining ones — positive. To minimize the value for energy 
one needs values of probabilities, corresponding to negative values of epsilons, as maximal 
as possible. Namely the possibility to obtain this formula for the eigenvalue of the many- 
particle Hamiltonian is basic to the Hall- Post lower bounds for energies construction [5|. 
This lower bound can be obtained taking a maximal value of probability corresponding to 
the lowest value of RH energy. This means: 

E Hall-Post = ^i- (27) 
Obviously, this is namely a lower bound because 

n n n 

E - Euau-Post = SiUi - £i^^i = ^ {^i - ^l) > 0. (28) 

i=l 1=1 i=2 
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On the other hand, our consideration has shown that this lower bound value for energy 
corresponds to wave function equal to the first basic component. This component is not 
antisymmetric, hence the corresponding eigenvalue is significantly lower than that for the 
antisymmetrized function. 

The real problem in our task is a huge amount of basic states while describing strongly 
correlated systems and as consequence — a very big rank of matrix of the antisymmetrizer, 
hence a large order of the Hamiltonian matrix. One needs the best possible values for 
energy and the high quality corresponding wave functions at the lowest possible price. The 
recommendations for successful realization of this program are as follows. 

First of all, one does not need the complete calculation of matrix of the antisymmetrizer 
A. As one may conclude from this matrix presentation, Eq.( l20i) . for complete matrix and 
its spectral decomposition construction only r linearly independent columns or rows of this 
matrix, i.e., matrices Z and Q, are necessary. These define the matrix of spectral decomposi- 
tion F, Eq.l pTl) . necessary for further calculations. During calculations of this part of matrix 
A one can control how far this calculation is from finish because due to the mentioned above 
condition the sum of squares of entries of column must equal the corresponding diagonal 
entry. When this condition is fulfilled with necessary precision, one can stop calculations 
for this column and start calculating the next column. The corresponding submatrix Z will 
show how many columns one needs. As it is required, this matrix must be nondegenerate, 
hence its determinant must be nonzero. Only in this case, one can continue column cal- 
culations. If the determinant equals zero, the basic state, corresponding to column under 
consideration, is linearly dependent on already calculated columns. This statement follows 
from relation for entries of projection matrix, stating that each submatrix of the projection 
matrix equals the Sturm matrix of corresponding rows or columns. Hence, the determinant 
of submatrix Z is the best indicator of the linear dependence of rows (columns) of projection 
matrix. Finally, one can check the rank of matrix. It must equal the trace. 

The last operation, when both matrices, F and H, are calculated, is construction of the 
Hamiltonian matrix and its diagonalization. This is the standard operation. However, in 
many applications only a few lowest eigenvalues and eigenvectors of this matrix are necessary. 
The method of obtaining converging approximations to exact values for a few lower states 
avoiding calculation of the complete Hamiltonian matrix is developed. The idea is based on 
applied transparent presentation for the Hamiltonian. As usual, the diagonalization of this 
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matrix requires the set of orthogonal transformations 



R+F+HFR = (FR)+ HFR. (29) 

In our approach the transformation for the entire matrix is not necessary - one can transform 
the matrix F in a necessary way. The needed transformations are products of (r — 1) 



Householder reflections 



101 ]. i.e., R = RiR2...Rr-i transforming the first, second and all 
the following rows of matrix F to convert it into the lower triangle form. To some extent 
this way of matrix transforming looks like the well-known QR decomposition of matrix. If 
the first row of matrix F is f^"*" = (/ii/i2---/ir) , the first Householder matrix, whose order 
equals r, is 

R, = 1-2!^. (30) 

The second Householder matrix is 

where fg*" = ( 722/23 •••/2r) is the upper (right) part of the second row of matrix F and the 
order of unit submatrix is shown. The last matrix of transformation equals: 



R-r-l — f if+ • (^2) 

Here = (fr-irfrr) ■ After these transformations are performed, the new Hamiltonian 
matrix F'+HF' appears, where F' = FRiR2...Rr--i. Obviously, the new matrix has the 
same eigenvalues as the initial one. Moreover, F''*'F' = 1 and F'F'+ = A. The matrix F' 
has zeros in the upper triangle and maximal possible diagonal entries, equal to norms of 
parts of entries of rows, involved in the Householder reflection, i.e., fj^fi, f2''f2,..., f^ifr-i, 
respectively. 

Now one can obtain that nonzero entry in the first row has only the first column of matrix 
F'. This entry squared equals the maximal weight (probability) of the lowest eigenvalue of 
RH. Any combinations of this column with remaining columns of matrix F' can only lower 
this probability. The nonzero weight (probability) of the first and second of the lowest 
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eigenvalues of RH can produce only the first and second columns, etc. The set of negative 
eigenvalues of RH, as mentioned, is finite. Let this number be equal k. Not more than {k — 1) 
Householder reflections are necessary to obtain an acceptable result because, as already 
mentioned, each column with zero entries, corresponding to negative RH eigenvalues, while 
combining with the first k columns, can only lower the probabilities of negative epsilons, 
hence a negative part of the eigenvalue for the total Hamiltonian. 

V. SUMMARY AND CONCLUSIONS 

Summarizing, let us present the list of consequent steps for the application of the proposed 
method. 

The first step is construction of the Reduced Hamiltonian operator given in Eq. (4) , and 
generation of the complete system of square-integrable functions for this operator. The 
spectrum of this Sturm-Liouville problem can be discretized by applying simple boundary 
conditions for eigenf unctions. 

The next step is construction of the system of basic functions for components, taking 
arbitrary antisymmetric functions for "spectator" and eigenfunctions of RH - for "actor" 
basis. The components must be arranged in the order corresponding to nonvanishing diag- 
onal entries of matrix H, i.e., eigenvalues of the RH operator. 

Now one can start calculation of the antisymmetrizer matrix using factorization, given 
in Eq.(15). This presentation is optimal for component basis because due to antisymmetry 
of "spectator" and "actor" functions both antisymmetrizers present on the left-hand side as 
well as both antisymmetrizers, present on the right-hand side, may be omitted. Obviously, 
one must calculate only matrices of two nontrivial operators - Pn-2,n and Pn-3,n-iPn-2,n- 
Later on, combining them with the identity matrix according to Eq.(16), one can construct 
the matrix of the antisymmetrizer. Due to simple structure of this matrix, thoroughly con- 
sidered in Sec. 3, one needs only a part of this matrix, necessary for matrix F construction 
according to formulas, given in Eqs. (20) and (21). By the way, one can start calculations 
of eigenvalues for the total Hamiltonian at any number of columns (rows) of the antisym- 
metrizer matrix, once it is not less than the number k of negative eigenvalues of the RH 
operator. The complete calculation requires r columns of matrix A, but if one stops at some 
other number of columns calculated, this is more or less successful approximation for the 
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rank of this matrix. 

The last step is appUcation of the first, second, third and following Householder reflections 
and diagonalizations of corresponding matrices, whose orders equal the number of reflection. 
Due to the variational character of these approximations, the results obtained, as usual, will 
converge to precise values from above. 

The application of the method for noninteracting particles is the simplest and most 
transparent. In this situation, RH is proportional to the Hamiltonian, characteristic of 
any Shell Model picture. It is the Hamiltonian for one particle, moving in the external 
fleld, multiplied by the number of particles N, hence its eigenvalues are well-known one- 
particle energies, multiplied by the same number. As it was mentioned, the matrix of the 
antisymmetrizer now is independent of any dynamics, and it splits down into a direct sum of 
projection matrices for different conflgurations. One is free to take spectral decompositions 
for every of these matrices separately. If particles are without any correlations, i.e., if 
the system can be described by one determinant, probabilities of all one-particle states in 
given conflguration equal hence the sum has a well-known form given in Eq.(26). In 

case when one must ensure good quantum numbers for such a system, the probabilities of 
different states of RH will equal the squared coefficients of fractional parentage, obtained 
after spectral decomposition of the antisymmetrizer matrix. 

Let us take for illustration of proposed ideas a simple model for ^'^F terms in conflg- 
uration of four identical fermions with spin equal 1/2. The order n of this matrix equals 8, 
and the rank r equals 3. The matrix F of spectral decomposition is: 

^ y^SjU ^/5/U -^3/560 -^33/140^ 

F+ = -^/l/l5 ^2/105 y/3/8 ^/l/56 ^JIJ^ -0760 ^3/28 ^11/84 

^ -^4/15 -v/5742 v^l/14 -v^l/lS -^5/48 -^rjjlVl sJlljM j 

(33) 

Here, to make the text easy, one presents the transposed matrix. Let diagonal entries £j 
of matrix of RH be: 

- 7; -4; 0; 0.1; 0.2; 0.3; 0.4; 0.5. (34) 

These eigenvalues are not far from the real situation. Existence of a few bound states of 
RH and all remaining states situated at positive energies, i.e., in the continuous spectrum, is 
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characteristic of strong interaction. As mentioned, in our approach the continuous spectrum 
is discretized due to boundary conditions chosen. 

The first state (row) of (the first column of F) matrix produces a very pure value 
for energy equal —0.645, the second - even worse, equal —0.374 due to the smaller total 
weight of basic states, corresponding to the first two negative eigenvalues of RH. The result, 
produced while applying the third state (row) of F+, is the best one. It equals —2. 129. 
This value is closest to the lowest eigenvalue of the total Hamiltonian matrix F+HF, which 
equals —2. 532. 

The result of the first Householder refiection, the first row of F'^ = RfF^, gives the 
eigenvalue equal to —2. 380. After the second Householder refiection applied, the lowest 
eigenvalue of the second order matrix equals —2. 532, hence at a given precision it coincides 
with the exact eigenvalue of the total matrix of the Hamiltonian. 

The above mentioned Hall-Post lower bound for this eigenvalue equals —7.000. The diag- 
onalization of the corresponding third order matrix in uncorrelated harmonic oscillator basis 
due to the core present in realistic nucleon-nucleon potential can give only positive values 
for the RH operator, hence a positive value for energy of the whole system. Obviously, our 
result is a real achievement for very complex interactions, such as recent realistic nucleon - 
nucleon potentials. 

Finally, let us draw some parallel between the developed approach and the well-known 
Hartree-Fock Self-Consistent Field method. In both cases first of all the building blocks 
for the system wave function construction are necessary. The common feature in both ap- 
proaches is construction of these blocks taking into account as many dynamic correlations 
as possible. In SCF these, due to the chosen scheme of simple antisymmetrization in given 
configuration, are the best possible one-particle functions. In our approach, when one can 
choose the refined way of antisymmetrization, these are the eigenf unctions of the reduced 
Hamiltonian operator. Both approaches at this stage have some attractive points. The SCF 
method applies determinants or CFP. The price for this simplicity is well-known unavoidable 
appearance of the phenomenological central field, the conversion of the linear problem to 
nonlinear equation for one-particle wave functions, and the essential impossibility to take 
into account some dynamic correlations while applying this method. In our method, all dy- 
namic correlations are taken into account in advance by using RH eigenf unctions, however 
now one can perform antisymmetrization only in infinite space. The next step in SCF is 
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iterations for the self-consistent field and optimization of corresponding one-particle func- 
tions. In our approach, this corresponds to the application of k Householder reflections and 
diagonalizations of the Hamiltonian matrix of the growing order. This procedure is equiva- 
lent to the lowest configuration definition in the SCF method. Finally, if SCF approach in 
some cases does not work, superpositions of configurations are necessary. This essentially 
complicates the consideration. In our case in analogous situation, one needs to take more 
columns of the antisymmetrizer matrix, no more serious problems, no new complex opera- 
tions. Our approach can be applied to the translationally invariant functions, while for the 
Self- Consistent Field approach such a modification is impossible. 

As a more realistic application let us comment on the modifications of traditional cal- 
culation in the harmonic oscillator basis for the three nucleon system pJJ]. The number 
of basic states (components) with the number of oscillator quanta up to E=44 is 16744. 
This equals the dimension of both - RH and antisymmetrizer matrices. The rank of the 
antisymmetrizer matrix, equal to the dimension of the total Hamiltonian matrix, is 5553. 
The ground state energy of triton in this approximation for the realistic nucleon-nucleon 



potential Reid93 [12l] equals -7.59 MeV (the exact value of ground state energy for this 
potential is -7.63). Among 5553 basic components only 131 have "actor" wave-functions, 
corresponding to negative eigenvalues of RH in two nucleon channels 1"*"0 (^5*1 — ^ Di) and 
0^1 (^5*0). Therefore, applying 130 Householder reflections for columns of matrix F before 
the total Hamiltonian matrix calculation one can transform it to lower triangle form with 
only 131 columns, having nonzero entries, corresponding to negative eigenvalues of the RH 
operator. Namely these entries produce nonzero probabilities of corresponding states. As a 
result, only these 131 basic antisymmetric functions (transformed columns of matrix F) are 
necessary for construction of the total Hamiltonian matrix, ensuring the mentioned result for 
the triton ground state energy. Obviously, transformation of dimension of the Hamiltonian 
matrix from 5553 to 131 is a real achievement of the proposed method, giving opportunity 
to simplify the realistic nuclear calculations. 
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